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Abstract 

DELPHI results are presented on the inclusive production of two {KKtt)^ states 
in the mass region 1.2-1.6 GeV/c^ in hadronic Z decays at LEP I. The measured 
masses (widths) are 1274±6 MeV/c^ (29±12 MeV/c^) and_1426±6 MeV/c^ (51±14 
MeV/c^) respectively. A partial-wave analysis of the (KKir)^ system shows that 
the first peak is consistent with the I^{J^^) = 0+(l++)/(0-+) ao(980)7r and 
the second with the I^{J^'^) = 0+(l++) fr*(892)ir + c.c. assignments. The 
total hadronic production rates per hadronic Z decay are (0.165 ± 0.051) and 
(0.056 ± 0.012) respectively. These measurements are consistent with the two 
states being the /i(1285) and /i(1420) mesons. 
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1 Introduction 

The inclusive production of mesons has been a subject of long-standing study at LEPI, 
as it provides insight into the nature of fragmentation of quarks and gluons into hadrons. 
So far studies have been done on the S'-wave mesons (both ^5*0 and ^Si) such as vr and 
p, as well as certain P-wave mesons /2(1270), i^r2(1430) and /2(1525) (i.e. ^P2) and 
/o(980) and ao(980) (i.e. ^Pq) [1-3]. Very little is known about the production of mesons 
belonging to other P-wave multiplets (i.e. ^Pi and ^Pi). For the first time, we present 
in this paper a study of the inclusive production of two J^*" = l"*""*" mesons, the /i(1285) 
and the /i(1420) (i.e. ^Pi). 

There are at least four known nonstrange isoscalar mesons [3], I^{J^^) = 0"'"(1+"'") 
and l'^{J^^) — 0+(0~''"), in the mass region between 1.2 and 1.6 GeV/c^, which couple to 
the decay channel {KKnf. These are the /i(1285), ?7(1295), /i(1420) and r/(1440). All 
are seen prominently in the peripheral production from Ti~p interactions [3], indicating 
that, despite their decay into (i^T^Tr)", they are mostly nn states, where n = {u,d}. 
There exist possibly two additional states, I^{J^^) =?-(l+-) /ti(1380) and I^{J^^) = 
0"'"(1"'"+) /i(1510), which may harbour a large ss content, as they are produced with 
considerable cross-sections in the peripheral reactions involving K~p interactions [3]. 
Given this complexity in the (KK-k)^ systems, it is important to find which resonances 
among these are readily excited in inclusive hadronic Z decays. 

The DELPHI data for this study are based on the neutral KKn channel in the reaction 

Z ^{K^K^7r^)+X'^ (1) 

Section 2 is devoted to the selection process for the event sample collected for this analysis. 
The KKtt mass spectra are studied in Section 3. It is shown that the selection of the 
events with low M{Kg K"^) mass is the crucial criterion to reveal the presence of two 
signals in the /i(1285) and /i(1420) mass regions. A partial- wave analysis, carried out 
to explore the spin-parity content of the two signals, is described in Section 4. The 
measurement of the production rates and differential cross-sections is presented in Section 
5. Conclusions are given in Section 6. 

2 Experimental Procedure 

The analysis presented here is based on a data sample of 3.4 million hadronic Z decays 
collected from 1992 to 1995 with the DELPHI detector at LEP. A detailed description of 
the DELPHI detector and its performance can be found elsewhere [4,5]. 

The charged particle tracks have been measured in the 1.2 T magnetic field by a set of 
tracking detectors. The average momentum resolution for charged particles in hadronic 
final states, A{l/p), is usually between 0.001 and 0.01, depending on their momentum as 
well as on which detectors are included in the track fit. 

A charged particle has been accepted in this analysis if its momentum p is greater 
than 100 MeV/c, its momentum error Ap/p is less than 1 and its impact parameter with 
respect to the nominal crossing point is within 4 cm in the transverse (xy) plane and 4 
cm/sin^ along the beam direction (2;-axis), 9 being the polar angle of the track. 

Hadronic events are then selected by requiring at least 5 charged particles, 3 GeV as 
minimum energy of the charged particles in each hemisphere of the event (defined with 
respect to the beam direction) and total energy of the charged particles of at least 12% 
of the centre-of-mass energy. The contamination from events due to beam-gas scattering 
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and to 7-7 interactions is estimated to be less than 0.1% and the background from r+ r~ 
events is less than 0.2% of the accepted events. 

After the event selection, in order to ensure a better signal-to-background ratio for 
the resonances in the K^K^n^ invariant mass system, tighter requirements have been 
imposed on the track impact parameters with respect to the nominal crossing point, i.e. 
they have to be within 0.2 cm in the transverse plane and 0.4 cm/sin 6' along the beam 
direction. 

identification has been provided by the Barrel Ring Imaging Cherenkov (BRICH) 
detector for particles with momenta above 700 MeV/c, while the ionization loss measured 
in the Time Projection Chamber (TPC) has been used for momenta above 100 MeV/c. 
The corresponding identification tags are based on the combined probabilities, derived 
from the average Cherenkov angle and the number of observed photons in the BRICH 
detector, as well as the measured dE/dx in the TPC. Cuts on the tags have been applied 
to achieve the best signal-to-background ratio, while rejecting e^, /x^, p and p tracks. A 
more detailed description of the identification tags can be found in Ref . [6] . In the present 
case, the K"^ identification efficiency (typically 50% over the kaon momentum range of 
this analysis [6]) has been estimated by comparing the (1020) to K^K" signal in the 
experimental data with a sample of simulated events generated with JETSET [7] tuned 
with the DELPHI parameters [8] and passed through the detector simulation program 
DELSIM [5] . Agreement within ±4% is observed between the data and the simulation. 

The Kg candidates are detected by their decay in fiight into tt+tt". The details of 
the reconstruction method and the various cuts applied arc described in Ref. [9]. Our 
selection process consists of taking the V^^s passing the standard criteria for quality of 
the reconstruction plus a mass cut given by 0.45 < M^n^n^) < 0.55 GeX/c^. 

After all the above cuts, only events with at least one KgK'^7r~ or KgK~7r'^ combi- 
nation have been kept in the present analysis, corresponding to a sample of 547k events. 



3 KgK'^Tr"^ Mass Spectra 

The KgK^TT^ invariant mass distribution is shown in Fig. la). Also shown in the 
figure is the same mass spectrum with a K*{892) selection (0.822 < M{Kti) < 0.962 
Ge\/c^), which would be appropriate if the decay of a resonance had proceeded through a 
ir*(892) intermediate state. Neither histogram shows a visible enhancement in the mass 
region between 1.2 to 1.6 GeX/c^. This is due to the enormous background in this mass 
region coming from the high number of K^K^^t:^ combinations per event (11 on average) 
in inclusive Z decays. The key to a successful study of the /i(1285) and /i(1420) under 
the circumstances is to select events with low [Kg K^) mass (Fig. lb)). This has the effect 
of selecting both the possible ao(980)^7r^ decay mode and, in case of A'*(892)^ + c.c. 
decay, the interference region of the two i^*(892) bands on the decay Dalitz plot, while 
reducing substantially the general background for the KKn system. Varying this cut 
on the Monte Carlo generated events suggests a mass cut M{Kg K^) < 1.04 Ge\/c^ to 
maximize both /i(1285) and /i(1420) signals over background. The application of this 
cut on the experimental data is shown in Fig. 2, where two clear peaks are now seen in 
these mass regions where the mass resolution is 8 and 9 Me^/c^ respectively. Based on the 
Monte Carlo generated event sample, we have verified that neither signal was a refiection 
of resonances whose mass is in the 1.0 to 1.5 Ge\/c^ range, such as the (1020) and the 
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iiri(1270) mesons nor was faked by a possible misidentification of kaons or pions coming 
from the decay of these resonances. 



8000 '- a) 

7000 - 

6000 - 

5000 - 

4000 - 

3000 - 

2000 - 

1000 - 



DELPHI 







1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 



M(K Kk), GeY/c 



■^o 3000 - 



^ 2000 

a 
> 

^ 1500 



b) Mjt#Hy,<( 



DELPHI 



t 2500 I * '"^''^'AV.A., 



1000 



1 1.05 1.1 1.15 1.2 1.25 1.3 1.35 /.< 



M(K K), GeV/c 



Figure 1: K^t^"^ a) and b) invariant mass distributions from the Z decays 

with the DELPHI detector at LEP I. The histograms with sohd circles are for the full 
data sample, that one with open circles is for data with a 0.822 < M{Kt:) < 0.962 GeV/c^ 
- K* selection. 



To estimate the background under the signals, we have used the Monte Carlo generated 
event sample from which we have removed all mesons with a major decay mode into 
(KKtt)^ in the mass region 1.2 to 1.6 GeV/c^. The resulting (KKn)^ mass spectrum was 
fitted between 1.15 and 1.7 GeV/c^ with a background function 

MM) = (M - Mo)"i exp(a2 M + M^), (2) 

where M and Mq arc the efi^ective masses of the (KKtt)^ system and its threshold, 
respectively, and are the fitted parameters (cti = 0.7 ± 0.1, 02 = 5.8 ± 2.0, — 
—3.1 ± 0.6). Then we have fitted the experimental {KKtt)'^ spectrum from 1.19 to 1.7 
GeX/c^ with the background function fb{M) determined above, adding two 5'-wave Breit- 
Wigner forms 

^"^^^^ " (M - Mry + (r,/2)2 

We have not used the relativistic angular-momentum dependent Breit-Wigner form in 
(3) because such a form would require a complete knowledge of the quantum numbers 
of the resonances as well as the branching ratios for all possible decay modes and the 
corresponding orbital angular momenta. The limited statistics, the high background 
under the signals are additional reasons to adopt the simple Breit-Wigner form (3). The 
fitted masses and widths, and F^, given in Table 1, are compatible with the PDG 
[3] values for the /i(1285) and /i(1420) resonances. The systematic uncertainty on the 
masses has been estimated from the various background fits (described later) to be about 
1-2 Mey/c^ and has been added quadratically to the statistical error. It should be noted 
that the parameters of the first peak are not compatible with those of the 77(1295). 
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Table 1. Fitted parameters and numbers of events 



Mass (MeY/c^) 


Width (MeY/c^) 


Events 


1274 ± 6 

1426 ±6 


29 ±12 

51 ±14 


358 ± 93 (stat) ± 59 (sys) 
870 ± 128 (stat) ± 136 (sys) 



The numbers of events in Table 1 correspond to the fit shown in Fig. 2, where the 
widths of the two peaks have been fixed to the fitted values, while the background pa- 
rameters were left free. 




1.2 1.3 1.4 1.5 1.6 1.7 



M(KK7t), GeV/c' 

Figure 2: Invariant mass distributions for the system tt^ with a mass cut 

M{Kg K^) < 1.04 Geyc^. The two sohd curves in the upper part of the histogram 
describe Breit-Wigner fits over a smooth background (sec text). The lower histogram 
and the solid curve give the same fits with the background subtracted and amplified by 
a factor of two. 

The main sources of systematic uncertainty come from the various cuts and selection 
criteria applied for the reconstruction and the charged K identification on the one 
hand and the conditions of the fit procedure on the other. To estimate the first type of 
error, we have compared the K^K^ mass distributions of the simulated sample with the 
real data. Normalized to the same number of events, the distributions agree within ±7%, 
in the low K^K^ mass region. 
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The /i(1285) and /i(1420) signals show up over a large background (~ 80%). Vari- 
ations of the background shape and amplitude induce sizable fluctuations of the fitted 
numbers of signal events. To quantify this effect, we have performed various series of fits, 
one varying the mass range of the fit, another leaving free the background parameters 
while fixing the width of the signals, another with a polynomial shape for the background, 
thereby allowing the background level and shape to fluctuate. In this way we estimate 
the uncertainty of the number of fitted events to be ±15% for the /i(1285) and ±14% for 
the /i(1420). The systematic uncertainties have been added quadratically and are shown 
in Table 1. 

The overall efficiencies for the two states have been estimated from the Monte Carlo 
simulated events to be: 

(0.063 ± 0.003)% for /i(1285), 
(0.45 ± 0.02)% for /i(1420). (4) 

The quoted numbers include the following corrections for the /i(1285) and /i(1420) 
respectively: branching ratios to KKn (0.09,1.), fractions of final states with charged 
pion (1/2,2/3), branching ratio of — > tt~^tt~ = 1/2 x 0.686 = 0.343, reconstruction 
and identification efficiency for the selected events (0.058, 0.061) and correction factor 
(0.70, 0.32) for the M{Kg K^) < 1.04 GeV/c^ mass cut. The quoted errors are statistical 
errors on the Monte Carlo sample. 

4 Partial- Wave Analysis 

In an attempt to get more information on the spin-parity content of the two signals we 
have performed a mass-dependent partial- wave analysis (PWA) of the K^K^t:'^ system. 
There have been many 3-body partial-wave analyses; the reader may consult PDG [3] for 
earlier references, for example, on ai(1260), 02(1320), fsTi (1270/1400) or i^2(1770). For 
the first time, we apply the same technique to a study of the {KKt:)^ system from the 
inclusive decay of the Z at LEP. 

A spin-parity analysis of the system composed of three pseudoscalars requires five 
variables, which may be chosen to be the three Euler angles defining the orientation of 
the 3-body system in its suitably-chosen rest frame and two effective masses describing 
the decay Dalitz plot. We have chosen to employ the so-called Dalitz plot analysis, 
integrating over the three Euler angles. This entails an essential simplification in the 
number of parameters required in the analysis, as the decay amplitudes involving the 
D-functions defined over the three Euler angles and their appropriate decay-coupling 
constants, are orthogonal for different spins and parities [10]. In these conditions, the 
mass-dependent PWA comes to fitting Dalitz plots, thus providing the contribution of the 
various 

jPC 

waves as a function of the M{KK7t). The actual fitting of the data is done by 
using the maximum-likelihood method, in which the normalization integrals are evaluated 
with the Monte Carlo events [11], thus taking into account the M{K^ K"^) < 1.04 G(N/c? 
cut. The comparison between fits is made on the basis of their maximum likelihood values 
and their description of the (KKn) , (Kn) and (KK) mass distributions. 

The first step of the analysis is to parametrize the background under the two signals. 
This background accounts for different processes with, for example, different overall multi- 
plicities. From a study of the side bands (away from the resonance regions), it has been de- 
termined that the background contains substantial amounts of ao(980) and K*{892) unas- 
sociated with the resonances. We have thus assumed that the background, which should 
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not interfere with the signals, can be described by a non-interfering superposition of a 
constant three-body phase-space term and the partial waves I^{J^^) = 0"'"(0^+) ao(980)7r 
(,5-wave)^ 0+(l++) (X*(892)X + C.C.) (,S-wave) and 0-(l+-) (X*(892)X + c.c.) (,S-wave). 
The {KKtt) mass dependence of the background components is parametrized by the 
phase-space-like form given by (2), but with = 0. 

The signals themselves are represented by a set of quasi two-body amplitudes which 
have the form of Breit-Wigner functions multiphed by spin-parity terms [11]. Those 
include the I'^[J^'^) = 0+(l"'"+), 0^(0"'^) and 0~(1"'"~) partial waves, where the possible 
decay channels ao(980)7r and K*{892)K + c.c. are allowed to interfere within a given J^'^. 

The PWA was performed for M^KKtt) in the 1.18 — > 1.66 GeV/c^ mass range. The 
first series aimed at determining the background contributions. In this fit, the signals 
were assumed to be composed of the 7^(1295) and /i(1285) for the first peak and for 
the second of the ?7(1440), /ii(1380) and /i(1420) resonances which were parametrized 
as Breit-Wigner forms with masses and widths fixed to the PDG values [3]. The fit was 
checked to well reproduce the (Kn) and (KK) mass distributions outside the regions of 
the peaks. 

The following step consisted in fitting the spin-parity content of the two signals. For 
this the individual background contributions were fixed to their fitted values, only the 
overall background rate was left free. The J^"" amplitudes were introduced individually to 
probe the spin-parity content of each signal with the mass and width of the corresponding 
{KK7r)° resonance being fitted. 

The results of these fits are the following. The I'^{J^'-'') = 0^(1^^) ao(980)7r and 
0^(0^+) ao(980)7r waves account equally well for the 1.28 GeV/c^ region, with the same 
maximum likelihood value, i.e the first peak is equally likely to be the /i(1285) or the 
7^(1295). On the other hand, if the mass and the width of the resonances are fixed to 
their PDG values, this region is better fitted by the I^(J^'^) = 0"'"(1"'"+) ao(980)7r wave. 
This refiects the fact that the first peak position is closer to the /i(1285) mass than to 
the 77(1295) mass, as already noticed in the fit of the M{KKtt) spectrum. 

In the 1.4 GeV/c^ region, the maximum likelihood value is better for the I'^{J^^) = 
0+(l++) K*{892)K + c.c. wave over the I^{J^^) = 0+(0-+), the 0-(l+-) K*{892)K+c.c. 
and the 0^(1'*"^) ao(980)7r waves by about 4, 8 and 9 units, respectively, and thus favours 
the /i(1420) K*K hypothesis over the 77(1440) or /ii(1380) ones, as is verified on 
the projections of the individual fits (not shown) on the {KKn), {Kn) and {KK) mass 
distributions. 

The results ofjhe best fit with the J^(J^^) = 0+(l++) ao(980)7r : /i(1285) and 
0^(1++) K*{892)K + c.c. : /i(1420) amplitudes are shown in Fig. 3 with the background 
contributions in the form of error-bands. The masses, the widths and the numbers of 
events found in the fit are statistically consistent with those given in Table 1. One ob- 
serves that the background events I^{J^^) = 0+(l++) {K* (892) K + c.c.) (S-wave) and 
0^(1+^) {K*{892)K + c.c.) (S-wave) exhibit a shape suggestive of a resonance. However 
this effect is simply the result of: i) our mass dependent function which decreases rapidly 
at high (KKn) mass to reproduce the fall off due to the M{Kg K^) < 1.04 GeV^c^ cut; 
ii) the K*{892)K threshold around 1.4 GeV/c^ whose effect is gradual because of the finite 
width of the K*{892). 

The major sources of systematic uncertainties come from the background description 
and the conditions of the overall PWA fit. To estimate them we have carried out series 
of fits leaving free the background contributions, then fixing the mass and width of 
the (KKn)^ resonances to their PDG values [3]. This was repeated with a polynomial 
background in M{KK'k) in place of the phase-space-like one. In all these fits, the partial 
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waves were fitted concurrently, in the form of a non-interfering superposition, to estimate 
their relative contribution. All the fits confirm the observations made previously. 

Taking into account the systematic uncertainties computed from the various fits, the 
number of /i(1285) and of /i(1420) events for the best fit of Fig. 3 are 237 ± 60(stat) ± 
70(syst) and 711 ± lOO(stat) ± 75(syst), respectively, consistent within systematic uncer- 
tainties with the values of the fit to the M^KKti) spectrum described in Section 3. 

All fits confirm the dominance of the I^{J^^) = 0+(l++) K*(892)i? + c.c. wave in 
the 1.4 Ge\^c^ region. The largest contributions of the 77(1440) and /ii(1380), estimated 
from the highest fitted rates of the J^(J^^) = 0+(0-+) and K*{m2)K + c.c. 

waves, correspond to production rates per hadronic Z decay of 0. +0.007 and (0.017l2;oi5) 
respectively, assuming a i^'*(892)^ + c.c. branching ratio of 100% for these resonances. 

5 Production Rates and Differential Cross-sections 

Prom the histogram fit described in section 3., we have measured the production rate 
{n) per hadronic Z decay for /i(1285) and /i(1420). The results are 



taking a KKn branching ratio of (9.0 ± 0.4)% for the /i(1285) and 100% for the 
/i(1420) [3]. 



The total production rates, per spin state and isospin, for the scalar, vector and tensor 
mesons with different strangeness, as a function of the mass [12,13] are shown in Fig. 4 for 
the averaged LEP data. To this figure we have added our measurements for comparison. 
It is seen that both /i(1285) and /i(1420) come close to the line corresponding to mesons 
whose constituents are thought to be of the type nn. This suggests that both /i(1285) 
and /i(1420) have little ss content. Indeed, the two states which arc thought to be pure 
ss mesons, the and the /2(1525), are down by a factor f« 1/4 (7 = 0.50 ± 0.02 [12] 
and k = 2, where k is the number of s and s quarks in the meson), as shown in Fig. 4. A 
high strange quark content is highly unlikely given the production rate (5). 

For completeness, we give in Fig. 5 and in Table 2 the /i(1285) and /i(1420) differential 
rates and cross-sections as a function of the scaled momentum Xp {xp = Pkk-kO /Phcam) for 
Xp > .05, as the signal to background ratio is too small for lower momenta. Quantitative 
comparison with JETSET predictions is not possible in a meaningful way as there was 
no tuning for /i(1285) and /i(1420) and the implementation of the {KKny decay of 
both resonances in JETSET had been done according to phase-space and not according 
to the correct spin-parity matrix element. The small excess of events in Fig. 5, near 
M{K^K^n^) = 1.55 GeV/c^ for 0.05 < Xp < 0.1, is not significant enough for further 
study here. 



in) 
(n) 



0.165 ±0.051 for /i(1285), 
0.056 ±0.012 for /i(1420). 



(5) 



Table 2. Measured production rates per hadronic event, differential cross- 
sections and experimental efficiencies for the /i(1285) and /i(1420), as func- 
tions of Xp. 



Xp interval 


/i(1285) rate 


{l/ah){da/dxp) 


Efficiency 


.05-. 10 
.10-.20 

.20-1.0 


0.046±0.026 
0.053±0.024 

0.051 ±0.022 


0.92±0.52 
0.53±0.24 

0.06±0.03 


(6.5 ± 0.7) X 10-4 
(9.4 ±0.8) X 10-4 

(6. 1 ± 0.7) X 10--^ 




Xp interval 


/i(1420) rate 


{l/ah){da/dxp) 


Efficiency 


.05-. 10 
.10-. 20 
.20-1.0 


0.018±0.006 

0.017±0.004 
0.015±0.005 


0.36±0.12 

0.17±0.04 
0.02±0.01 


(3.1 ±0.3) X 10-=^ 

(8.5 ±0.5) X 10-3 
(3.9 ±0.3) X 10-3 
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Figure 3: M^K^ i^^vr^) distributions per 10 Me\^c^ with a breakdown into the partial 
waves for the signals in the lower histogram and for the background shown as one er- 
ror-band. The signals consist of 1++ ao(980)7r for the first peak and 1++ ii'*(892)^ for 
the second peak. The background is composed of a non- interfering superposition of (1) 
isotropic phase-space distribution and the following partial waves: (2) 0~+ ao(980)7r, (3) 
1++ K*{SQ2)K and (4) l+~ ii'*(892)^ waves. The last three contributions are shown 
magnified by a factor of two. 
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Figure 4: Total production rate per spin state and isospin for scalar, vector and ten- 
sor mesons as a function of the mass (open symbols) from [12]. The two solid circles 
correspond to the /i(1285) and the /i(1420) measurements presented here. 
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Figure 5: The tt^ invariant mass spectra for various Xp intervals as indicated. Dots 

are the data, the sohd hues show the result from the fit and the background contribution. 



12 



6 Conclusions 

We have studied the inclusive production of two (KKn)^ states in Z decays at LEP I. 
The measured masses and widths are 1274 ± 6 and 29 ± 12 MeV/c^ for the first peak and 
1426 ± 6 and 51 ± 14 MeV/c^ for the second one, compatible with those of the /i(1285) 
and /i(1420) mesons [3]. For the first time, a partial-wave analysis has been carried 
out on the (KKn)'^ system from the inclusive Z decay While the results are ambiguous 
between the I^'{J^^) = 0+(l++) and 0+(0-+) ao(980)7r waves in the 1.28 GeV/c^ region, 
the second peak is uniquely consistent with I^{J^^) = 0"'"(1"'""'") i^'* (892)^ + c.c. On the 
other hand, the comparison of the hadronic production rate of these two states with a 
previous study of the production rate [12,13] for the S — 1 mesons (which included ^Si, 
^Pq and ^P2) suggests that their quantum numbers are very probably I^{J^'-') = 0+(l++) 
and that their quark constituents are mainly of the type ran, where n = {u, d} and thus 
confirms that these states are very likely the /i(1285) and /i(1420) mesons. Finally, we 
conclude that the mesons ?7(1295), ?7(1440) and /ii(1380) are less likely to be produced in 
the inclusive Z decays compared to the /i(1285) and /i(1420). 

Acknowledgements 

We are greatly indebted to our technical collaborators, to the members of the CERN- 
SL Division for the excellent performance of the LEP collider, and to the funding agencies 
for their support in building and operating the DELPHI detector. 
We acknowledge in particular the support of 

Austrian Federal Ministry of Education, Science and Culture, GZ 616.364/2-III/2a/98, 
FNRS-FWO, Flanders Institute to encourage scientific and technological research in the 
industry (IWT), Federal Office for Scientific, Technical and Cultural affairs (OSTC), 
Belgium, 

FINEP, CNPq, CAPES, FUJB and FAPERJ, Brazil, 
Czech Ministry of Industry and Trade, GA CR 202/99/1362, 
Commission of the European Communities (DC XII), 
Direction des Sciences de la Matiere, CEA, Prance, 

Bundesministerium fur Bildung, Wissenschaft, Forschung und Technologic, Germany, 
General Secretariat for Research and Technology, Greece, 

National Science Foundation (NWO) and Foundation for Research on Matter (FOM), 
The Netherlands, 
Norwegian Research Council, 

State Committee for Scientific Research, Poland, SPUB-M/CERN/PO3/DZ296/2000, 

SPUB-M/CERN/PO3/DZ297/2000 and 2P03B 104 19 and 2P03B 69 23(2002-2004) 

JNICT- Junta Nacional de Investigagao Cientifica e Tecnologica, Portugal, 

Vedecka grantova agentura MS SR, Slovakia, Nr. 95/5195/134, 

Ministry of Science and Technology of the Republic of Slovenia, 

CICYT, Spain, AEN99-0950 and AEN99-0761, 

The Swedish Natural Science Research Council, 

Particle Physics and Astronomy Research Council, UK, 

Department of Energy, USA, DE-FG02-01ER41155, 

EEC RTN contract HPRN-CT-00292-2002. 

S.-U. Chung is grateful for the warm hospitality extended to him by the CERN staff 
during his sabbatical year in the EP Division. 



13 



References 

[1] DELPHI CoUab., P. Abreu et ai, Phys. Lett. B449 (1999) 364. 

[2] OPAL Collab., R. Akers et ai, Z. Phys. C68 (1995) 1. 

K. Ackcrstaff et ai, Eur Phys. J C4 (1998) 19; ibid. C5 (1998) 411. 

[3] Particle Data Group, K.Hagiwara et al., Phys. Rev. D66, 010001 (2002). 

[4] DELPHI Collab., P. Aarnio et al, Nucl. Inst. Meth. A303 (1991) 233. 

[5] DELPHI Collab., P. Abreu et al, Nucl. Inst. Meth. A378 (1996) 57. 

[6] DELPHI Collab., P. Abreu et al, Eur. Phys. J. C5, (1998) 585. 

[7] T. Sjostrand, Comput. Phys. Comm. 82 (1994) 74. 

[8] DELPHI Collab., P. Abreu et al, Z. Phys. C73 (1996) 11. 

[9] DELPHI Collab., P. Abreu et ai, Z. Phys. C65 (1995) 587. 
[10] See, for example, S. U. Chung, 'Spin Formalisms,' CERN preprint 71-8 (1971). 
[11] E852 Collab., S. U. Chung et ai, Phys. Rev. D60 (1999) 092001. 
[12] V. Uvarov, Phys. Lett. B511 (2001) 136. 

V. Uvarov, Phys. Lett. B482 (2000) 10. 
[13] P.V. Chhapnikov, Phys. Lett. B525 (2002) 1. 



